Adiabatic nonlinear waves with trapped particles: II. Wave dispersion 
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A general nonlinear dispersion relation is derived in a nondifferential form for an adiabatic si- 
nusoidal Langmuir wave in coUisionless plasma, allowing for an arbitrary distribution of trapped 
electrons. The linear dielectric function is generalized, and the nonlinear kinetic frequency shift a;NL 
is found analytically as a function of the wave amplitude a. Smooth distributions yield wnl oc y^, as 
usual. However, beam-like distributions of trapped electrons result in different power laws, or even 
a logarithmic nonlinearity, which are derived as asymptotic limits of the same dispersion relation. 
Such beams are formed whenever the phase velocity changes, because the trapped distribution is 
in autoresonance and thus evolves differently from the passing distribution. Hence, even adiabatic 
i^NL(a) is generally nonlocal. 
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I. INTRODUCTION 

Within the geometrical-optics approximation, adia- 
batic waves in coUisionless plasma are described conve- 
niently within the average-Lagrangian formalism orig- 
inally proposed by Whitham [11 H]. In our Ref. [3], 
further called Paper I, this formalism was restated to 
also accommodate effects of particles trapped in wave 
troughs. Specifically, the wave Lagrangian density £ was 
derived explicitly under the assumption that the number 
of trapped particles within each wavelength is a fixed in- 
dependent parameter of the problem. Hence the wave 
nonlinear dispersion relation (NDR) can be expressed 
in a nondifferential form which is more easily tractable 
than solutions obtained straightforwardly from integrat- 
ing the Vlasov-Maxwell system (see, e.g., Refs. [IHS]) 
or those inferred from the "multiwater bag" formulation 
[lU] . Namely, the general NDR for an adiabatic wave 
reads as (Paper I) 



da^ - 0, 



(1) 



a being the wave amplitude; thus, analyzing the function 
das' would be sufficient to infer the wave dispersion. 

Here, we perform a systematic study of kinetic nonlin- 
earities that Eq. (IT]) yields for one-dimensional sinusoidal 
electrostatic waves with trapped particles, extending our 
brief letter [111 and providing important technical details 
missing in that work. Specifically, a general NDR is de- 
rived in a nondifferential form for an adiabatic sinusoidal 
Langmuir wave in coUisionless plasma, allowing for an ar- 
bitrary distribution of trapped electrons [Eq. (111]. The 
linear dielectric function is generalized, and the nonlinear 
kinetic frequency shift wnl is found analytically as a func- 
tion of the wave amplitude a. Smooth distributions yield 
wnl oc ^/a, as usual. However, beam-like distributions of 
trapped electrons result in different power laws, or even 
a logarithmic nonlinearity, which are derived as asymp- 
totic limits of the same dispersion relation. Such beams 
are formed whenever the phase velocity changes, because 
the trapped distribution is in autoresonance and thus 
evolves differently from the passing distribution; hence. 



even adiabatic WNL(a) is generally nonlocal. These re- 
sults are also used in Paper 111 [12 , where evolution of 
waves with trapped particles is studied. 

The work is organized as follows. In Sec.|TT] the general 
NDR is derived, connecting the wave frequency a; with 
the wave number k for a given a and particle distribu- 
tion; also, the properties of the NDR kernel are explored. 
In Sec. |III[ the general NDR is given an asymptotic rep- 
resentation for distributions smooth near the resonance; 
in particular, the linear dielectric function is generalized. 
In Sec. |IV[ the effect of particle beams on NDR is stud- 
ied. In Sec.|V] implications of these results are discussed; 
specifically, WNL(a) is predicted to be a nonlocal function, 
and the friction drag on trapped particles is predicted to 
affect the wave frequency sweeping in coUisional plasmas. 
In Sec. |V1[ we summarize our main results. Some auxil- 
iary formulas are also presented in the appendix. 



II. GENERAL DISPERSION RELATION 

A. Wave Lagrangian 

From Paper 1, it follows that the Lagrangian density for 
a one-dimensional adiabatic nonlinear electrostatic elec- 
tron mode in coUisionless plasma reads as 



£ = £e 



Stt 



.nW(H(^))/^ 



i(*)m(*) 



(^^'^)/., (2) 



assuming that the ion contribution is negligible. Here, 
£em describes quasistatic fields (if any), the second term 
is the time-averaged energy density of the wave electric 
field E, n^P'> and n'^*^ are the densities of passing and 
trapped electrons, respectively, and 'H'-''-' and "H*-*-* are the 
oscillation-center (OC) energies averaged over the corre- 
sponding distributions fp and ft of the particle canonical 
momenta. Specifically, 



H(p) = £ + Pu - mu^ /2, 
UW =£- mu^lX 



(3) 
(4) 



where m is the electron mass, and £( J, a) is the electron 
energy in the frame K moving with respect to the lab- 
oratory frame K at the wave phase velocity u — uj/k. 
Also, J is the trapped-particle canonical momentum, or 
the canonical action in K, and 



P ~ mu ± kj 



(5) 



is the passing-particle canonical momentum, where the 
plus sign is chosen for particles traveling with velocities 
larger than u and the minus sign is chosen otherwise. 
Combining Eqs. ([2|-(|4]), one gets 



£ = £n 



-nlP\P)f^u- 






-n{£)j 



where we introduced the electron total density n = n^P^ + 
rv-^', and the energy £ is now being averaged over all 
particles. By definition (Paper I), the particle densities, 
oj, and k are arguments of £ that are independent of a. [It 
is only through solutions of Eq. (IT]) that they will become 
interconnected.] Hence, Eq. (IT]) yields the general NDR 
in the form 



da{E^)/{8TT)-nda{£)f^Q, 



(6) 



where da is taken at constant (w, k) and also at fixed 
canonical momenta of individual particles. Specifically, 
fixed at differentiation can be J, not only for trapped 
particles but also for passing ones, due to Eq. (Is]). 



B. Sinusoidal- wave approximation 

Equation (|6| describes the dispersion of one- 
dimensional electron Langmuir waves completely, ac- 
counting simultaneously for linear and nonlinear disper- 
sive effects, both fluid and kinetic. In principle, each 
of those effects can then be assessed from the general 
NDR by substituting the appropriate £{J,a) [13]. Yet, 
we will address only one specific regime here, as an illus- 
tration of how our Lagrangian formalism can be useful 
for finding NDRs explicitly; namely, we will assume that 
the wave can be considered monochromatic. The condi- 
tions under which this assumption is justified as a precise 
asymptotic limit at small amplitudes are discussed in de- 
tail in Refs. [HITS] (see also Refs. [HUT]) and will not 
be repeated here. Instead, our purpose will be to show 
that, once the sinusoidal- wave approximation has been 
adopted (like in traditional calculations) , our Lagrangian 
formalism predicts new effects and also yields simulta- 
neously various types of kinetic nonlinearities that were 
previously known from different contexts. 

In a sinusoidal wave, the oscillations of both passing 
and trapped electrons in K satisfy 



p^ /{2m) + e(po cos{kx) = £, 



(7) 



where x and p are the particle coordinate and momen- 
tum, and e is the electron charge. (For clarity, we will 



(a) passing 




(b) trapped 




FIG. 1; Schematic of single-particle trajectories in phase 
space {x,p), illustrating the definition of 2nJ (shaded area): 
(a) for a passing particle, (b) for a trapped particle. For J to 
be continuous at the separatrix (dashed), with the passing- 
particle action defined as 2nJ = f pdx, for trapped particles 
one must use the definition 2ttJ = {1/2) § pdx. This area 
is encircled by a particle within half of the bounce period; 
thus, the corresponding canonical frequency Q, equals twice 
the bounce frequency. 



assume e(^o > and fc > 0. Also, the large-scale elec- 
trostatic potential (p is omitted here because it does not 
affect the wave local properties, albeit it may affect the 
dynamics of the wave (Paper III) and of the particle 
distribution [H]-) The canonical action J is hence in- 
troduced as follows. For passing particles we will use 
the definition J = {2Ti)^^§pdx, and thus J ~ \p\/k 
for large p. For trapped particles, however, we will use 
J = [An)'"^ § pdx, so that J is continuous across the 
separatrix (Fig. IT]) . From Eq. (J7| , one can write then 



J = ^/a jj. 



(8) 



Here we introduced the wave amplitude as a = 
eip^k'^ / {muj"^) > 0; also, J = muj/k'^, and the dimension- 
less action j is a continuous function of the normalized 
energy r ^ {£ + e(fo)/{2e(po), 



J = i{r)- 



(9) 



Specifically, j — for a particle resting at the bottom of 
the wave trough (r = 0) , with the corresponding value at 
the separatrix (r = 1) being j* = 4/7r (Appendix [A]) . 
Hence, one can express £ as 

£{J,a) = [2ar{j)~a]Juj, (10) 

with j{J,a) taken from Eq. (|8|. Then Eq. (|6| reads as 



UJ 



GF{J) dJ, 



(11) 



where 



A-Kne^ /m is the electron plasma frequency 



squared, F{J) is the action distribution normalized such 
that j^ F{J) dJ = 1, and 



G = 



da£{J,a) 



(12) 



Equation ( 11 ) is a master equation, which accounts for 



all dispersive effects simultaneously, both linear and non- 
linear (to the extent that the sinusoidal-wave approxima- 
tion applies |19|). In particular, it shows that the contri- 
bution to w^ of particles with given J is determined solely 
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FIG. 2: Weight fun ctio n G(j). Dashed are approximate solu- 
tions given by Eq. ( 19 1, the asymptote, and j = j*. 



by the kernel G, thus acting as the weight function. We 
study the properties of this function below. 



C. Weight function G 



From Eq. (12), one gets 



G = [da{2ar - a)]j = 2r-l + 2a (9,r)„ 



(13) 



where the partial derivative is taken at fixed J (as de- 
noted by the subindex). The latter term is found from 



{dar). 



{daJ)r 
{drJ)a 



J 

'2aj'' 



(14) 



where we substituted Eq. ([8|. Thus, G — 2r — 1 - j/f, 
where the right-hand side here is a function of r only, or 
a function of j only [due to Eq. ^]. We will introduce 
these two representations separately, with 

g{r) = 2r-l-j{r)/j'{r) (15) 

(Appendix ro serving as an auxiliary function for 

G{j) = 9{r{j)). (16) 

Hence, the properties of G{j) are understood as follows. 
First, let us consider asymptotics of G{j), starting with 
j ^ 1. As seen from Eq. ( [T5| , G(0) = — 1, which also 
flows from the fact that £ w ~e(pQ for particles close to 
the bottom of the wave trough. The next-order term can 
be inferred from £ « — e<^o + J^^ which is due to oscil- 
lations of deeply trapped particles being approximately 
linear. In particular, one can take fi « 2loe (Fig. Ilj, 
where uje — {ek'^cpo/m)^''^ = uJ^/a is the linear bounce 
frequency HO]- Hence, Eq. (12) yields G •■ 



-1+J. 



The asymptotic behavior of G{j) at j » 1 is un- 
derstood by taking £ « p^ /{2m) + $, where $ = 
e^TOi^Q/(4p^) is the ponderomotive energy in the dipole 
approximation (Paper I) for a zero-frequency wave 



[Eq. (It])]. Remember that \p\ w Jk in this case, so 
{da£)j = 2$/a. Thus, 



G(j » 1) « l/(2j-2). 



(17) 



and higher-order corrections can be found similarly. 

While negative at small j, G{j) is positive at large j. 
(This means that deeply trapped particles reduce w^, 
whereas those near the separatrix and untrapped ones 
increase w^; cf. Ref. [HI [55]).) Now let us consider its 
behavior near the separatrix, which is where the canoni- 
cal frequency i7 approaches zero. Since £ equals the OC 
energy in K (Paper I), one has 



rj = dj£ — Lo dj{2ar — a) = 2auj/j . 



(18) 



Therefore j'{r — 1) is infinite, and Eq. (15) yields that 
G be continuous (albeit nonanalytic), with G{j^) = 1. 
Hence, we summarize our results as follows: 



GU) = 



1, 

1 5 



2p 16je 



j » 1. 



(19) 



These are matched (Fig. [2]) by direct calculations in 
which an explicit expression for j(r) is substituted into 
Eq. (15). Also, manipulating that and related formulas 
given in Appendix [A] shows that the NDR arising from 
the waterbag formulation (Ref. [10], with zero driving 
field) can, in principle [33], be reduced to our Eq. (11). 



D. Small-amplitude cold-plasma limit 

It is instructive to consider, as an example, the cold- 
plasma case in the limit of small amplitudes. Then, all 
particles are passing, and J — \p\/k, so (Paper I) 



u ± kJ/m, 



(20) 



are 



where the sign index denotes sgn (w — u) . Yet v 
not present then, and v~ <C m, in which case Eq. (20) 
gives J •» J. [In other words, one may assume F{J) w 
5{J — J).] From Eq. (17), one thereby obtains G = a/2, 
and thus Eq. (11) predicts w^ = w^, as expected. 



Below we will calculate w^ also for other representative 
cases, by formally constructing various asymptotic ex- 
pansions of the integral in Eq. ( [TT] ). Remarkably, Eq. ( 19 ) 
will be sufficient for that; i.e., knowing the exact form of 
G(j) is not needed, except for evaluating insignificant 
numerical coefficients (Appendix [A|) . 



III. SMOOTH DISTRIBUTIONS 

First, let us consider the case when the distribution 
function F(J) is smooth, e.g., thermal. Of course, the 
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FIG. 3: (a) Function *(j), Eq. (221. (b) Function Q{j) 



Eq. (261. The dashed hnes show j = j, and also the asymp- 
tote of Q(j). For details see Appendix [A| 



wave dynamics could not be described then within the 
present approach, on the score of particle unavoidable 
trapping and detrapping that would be associated with 
variations of a, rendering the wave nonadiabatic. Yet, 
assuming that the wave is quasihomogeneous and qua- 
sistationary (which are our necessary requirements in 
any case), its local dispersion must be approximately the 
same as that of a truly homogeneous and stationary wave 
|24j . Since the latter has rS*''^ fixed and is clearly adia- 
batic, our formalism is hence applicable for deriving the 
local NDR and can be built on as follows. 

Suppose that F{J) remains finite at small J or, at 
least, diverges less rapidly than J^^. Then, one can take 
the integral in Eq. (|11[) by parts and obtain 



H 



T{J,a)F'{J)dJ^O, 



(21) 



where T(J, a) — '^{j)Jy/a, with the dimcnsionless func- 
tion \I/(j) [Fig.jSJJa)] introduced as 



*(j) 



G{j)d3. 



(22) 



At a <C 1, T changes rapidly with J compared to F{J). 
Then, without using the explicit form of T( J, a) but 
rather drawing on the leading terms in 



T(J,a) 




J^Ja^/^, 



. , J » Jai/2 



(23) 



[obtained from Eq. (19)], one can put Eq. (21) in an 



asymptotic form, specifically as follows. 



A. Asymptotic representation 

To approximate the integral 

2 



U 



T{J,a)F\J)dJ, 



let us replace TF' with 

TF' w iTF')i + {TF')2 



{TF') 



(24) 



(25) 



Here, the individual terms on the right-hand side corre- 
spond to asymptotics of TF' at small, large, and inter- 
mediate J, correspondingly, so that TF' is approximated 
accurately at all J. Specifically, we take 

(T{J)[F'{0) + F"{0)J], * = 1, 

(TF')»=| [ajV(2J)] F(J), J = 2, 

[ [aJV(2J)] [F'(0) + F"iO)j] , ^ - 3. 

[Due to T - aJV(2 J) = 0( J"'^) at large J, higher-order 
corrections are insignificant compared to those discussed 
below.] Hence, U = UiF"{0) + U2, where 



2 

Ui = - 
a 



^"|T(J,a)J-"^|dJ, 



U2 = lf {t( J, a)F'(O) - ^ [F'{J) - F'{0)] | dJ. 
It is convenient to rewrite Ui in terms of the function 



Q{j) = 1 - 2jvl/(j) 
Fig. ISFb)], which we will also use below. Namely, 



f/i 



-kJ \fa, K -- 



Q{z) dz. 



(26) 



(27) 



The value of the numerical constant >f can be readily esti- 
mated from FigJsFb) as xr ~ 0.5; for a precise calculation 
see Appendix I A| 

To simplify the expression for U2, notice that 



dU2 
da 



= 2F'(0) / 
Jo 



°° d_ 
da 



T(J,a) 



dJ, 



(28) 



which can be checked directly. On the other hand, 

"T(J,a) 



d_ 
da 



.P dj^j) dJiJ,a) 
J dj da 



J2 di^j) dj{J,a) 



2a dj 



and also 



dJ 



T{J,a)J 



2a2 



dJ 



dJ = 



d_ 
dJ 



T{J,a)J 
2a2 



T(J,a)J 



2a2 



4a' 



(29) 



(30) 



This yields 

daU2 = -J^F'(0)/{2a). (31) 

Integrating the latter with respect to a gives 

U2^ -[J^F'{0)/2]lna + Uc. (32) 

The integration constant Uc, independent of a, equals 
U2{a = 1), so it can be written as 



C/n 



2T{J,1)F'{0)+'^[F'{J)~F'{ 
I ■J 



Using that T( J, 1) = J*( J/ J), one gets 

'J 



Ur = J' 



F'{J)-F'{0)Q 



J 



dJ 



dJ. 



(33) 



[Notice that, although determined by essentially nonlin- 
ear dynamics in the narrow vicinity of the resonance, 
Q{J / J) nevertheless can affect the integrand on the ther- 
mal scale.] Finally, Eq. (24 1 is summarized as 



U = -[pF'(0)/2] In a - K^/aJ^F"[()) + U^- (34) 



B. Nonlinear frequency shift 



Using Eqs. (33) and (34), we can rewrite Eq. (21) as 

,2 



2fc2 






e(w,fc)H--^Cilna+-^xC2Va = 0. (35) 



Here we introduced 



e = 1 



m^ujp '■°° 



fc6 



F'(J)-F'(0)Q 



dJ 
T 



(36) 



and also 

Ci = (m/fc)2F'(0), C2 = im/k)^F"iO). (37) 

In particular, when F'(0) — 0, the nonlinear part of 
Eq. ( 35 ) is small, yielding that the nonlinear frequency 



shift ciJnl is also small; hence. 



-1 ,^, ,2 



de\ >€uj; /e(^o 
WNL = -l^| -nr\l C2. 



(38) 



Yet, at nonzero -F"(0), the nonlinear part of Eq. (35) 



diverges logarithmically at small a; i.e., wave interaction 
with resonant particles has a strong effect on uj. 



C. Comparison with existing theories 



Equations (35)- (38 1 generalize the existing NDR for 



quasimonochromatic eigenwaves in plasmas with smooth 
distributions of particles (THl [HI [HI 1^51 - 150] , namely, as 



follows. First of all, notice that e, serving as a general- 
ized linear dielectric function, is a functional of the action 
distribution. Unlike the commonly used distribution of 
"unperturbed" velocities fo{v), which depends on how 
the wave was excited [5T], F{J) is defined unambigu- 
ously; thus, the above equations hold for any excitation 
scenario. (Notice that finding F{J) itself is kept as a 
separate problem, which is also where the effect of the 
quasistatic field [H] comes in.) Second, even if put in 
terms of fo{v), Eqs. (35l-(38) cover a wider class of par- 
ticle distributions. The latter is seen as follows. 

For example, consider a homogeneous wave developed 
with fixed u slowly from zero amplitude, so each J is 
conserved [31], even through trapping and untrapping 
[H [MSU- Then F{J) = Fo{J), index henceforth 
denoting the initial state. Yet, since there was no wave 
in that state, Eq. (20) applies, so each ^th derivative of 
Fo(J) reads as 



rW, 



i\f(^)f.,+ 



Wi 



F^^'iJ) = {k/mY [f^^>{v+) + {-irf^'>{v-)\ . (39) 

Let us use bars to denote limits /g (w — >■ u±) and intro- 
duce 



fViv) 



fo (u-), v<u, 



(40) 



which is a piecewise-constant function equal to the left 
and right limits for v < u and v > u correspondingly; in 
particular. 






(41) 
(42) 



Hence e rewrites as follows: 



/;(«) - Q!«{v) 



where Q = Q{\v/u — 1|). (Remarkably, contributing to e 
are both passing and trapped particles.) 
First, compare Eq. (43) with the usual ^^l 



^--'-if' 



foiv) 



V — u 



dv, 



(44) 



P denoting the principal value [3S]. For smooth fo{v), 
our e can be put in the same form as e^j. because 



Qfl){v)dv = 0. However, Eq. (43) is valid 



P r {v~ 
J— 00 ^ 

also when /o(f) is discontinuous across the resonance, a 
case in which eL is undefined. This is because the inte- 



grand in Eq. (43) is finite (piecewise-continuous), so the 



integral converges absolutely rather than existing only 
as a principal value (like eL does). Second, for smooth 



/o(w), when Ci = and C2 = 2/^'(u), Eq. (g for wnl 
matches that in Ref. [21j, including the coefficient, since 
>c w 0.544 ( Appendix [A|. Yet, unlike the existing theory, 
our results apply just as well for arbitrary Ci and C2, in 
which case fo{v) may not be smooth while F{J) is. 



IV. BEAM-LIKE DISTRIBUTIONS 

Suppose now that, in addition to a smooth distribu- 
tion, near the resonance there is a phase-space clump or a 
hole, further termed uniformly as a beam with -?&( J) ^ 
and some average spatial density n^ ^ 0. Since particles 
contribute additively to the right-hand side in Eq. (Ill, 
wG immediately obtain 



2 2 

e(w, k) + —^ Ci Ina -t- -p- ycCi\J a 



^ = 0, 



where *8 is the beam-driven term. 



» 



2c.g 



anijUj^ Jq 



GFb{J)dJ, 



(45) 



(46) 



with w^ = Airnhe^ /m. [Unlike F{J), the distribution 
Fh{J) is normalized to the beam density here.] For ex- 
ample, consider a beam of deeply trapped particles, so 
that Fh{J) — nbS{J). Since G(0) = —1, one then has 



«8 = 2Lul/{auj^). 



(47) 



Hence, when the beam nonlinearity dominates, 

e{uj, k) + 2ojI/uI = 0. (48) 

For instance, within the model [37] 



e(a;,fc) = 1 



oj'^ — ik'^v^ ' 



(49) 



where vt is the thermal speed, one thereby gets the fol- 
lowing dispersion relation: 



w 



2„,2 



1 + 2^2/^1 



-f 3r wf. 



(50) 



in agreement with Ref. |15j . Strictly speaking, the 
sinusoidal- wave approximation required for this result to 
hold applies only at loI/uj'^ <C 1 [IS], so one can further 
approximate Eq. \hQ\ as 



a;2 ^ujI{1~ 2ul/uj%) + ik^v. 



Ti 



(51) 



in agreement with Ref. [38 . Since Eq. ( |49[ ) implies kvT <C 
w, one has lu fn ujp, so the above result also rewrites as 



2 2 o 2 / 



(52) 



where w^ 



jp + 3k^v^ corresponds to the linear limit. 
In particular, notice that wnl = 0{a^^), with cjnl < 
for a clump and wnl > for a hole. (Remember that 
w^ ex nf, may have either sign.) 

Now let ni, itself depend on the wave amplitude. For 
example, a Van Kampen mode would have w^/a in- 
dependent of a, yielding that trapped particles result 
in a linear frequency shift wnlj in agreement with the 
original linear theory [3S]- Also, let us consider the 



case when Fi, is constant across the trapping width, i.e., 
Fb{J) = F{,6(J)6(J* - J), _with Q being the Heaviside 
step function. Hence, rif, = Fi,J^,, and 



«B 



2^b' *(j* 



8FfcJ t^p 



STrny^ oj^ 



(53) 



where we used ^(j*) = 4/(37r) (Appendix [A]). Hence, the 
approximate formula is obtained, like above: 



U! 



^l 



SupFb 
37r-y/an 



(54) 



Equation (54) also matches the result found previously, 



e.g., in Ref. 'XL 



V. DISCUSSION 



The above analysis shows that the power index /3 in 
the scaling for the nonlinear frequency shift wnl oc a^ 
depends critically on the properties of F{J) at small J. 
The commonly anticipated /3 = 1/2 is realized when F{J) 
is smooth, an example being the case when the wave de- 
velops adiabatically from zero amplitude. However, per- 
turbations to F on the trapping-width scale J* result in 
additional nonlinearities which yield f3 < 0. In particu- 
lar, this affects the sign of da(^Nh, which can be decisive 
for the wave stability (Paper HI). Even more generally, 
describing such nonlinear kinetic effects in a fluid-like 
manner [3DJ [301 HI] is possible, in fact, only in special 
cases (e.g., when a{t) is monotonic and u is fixed, as in 
Ref. [33), because wnl is not a local function of a. 

For example, consider a driven nonlinear mode with 
adiabatic chirp imposed [IHl H5H45) . Due to the chirp, 
the wave phase velocity changes, carrying trapped parti- 
cles along because of the autoresonance; then F{J < J*) 
remains intact. Yet F{J > J*) is changing, because pass- 
ing particles do not conserve their actions as defined here 
^31,. Hence, even for an initially smooth F{J), the chirp 
generally produces discontinuity of the action distribu- 
tion at J = Jif, thereby changing the scaling for ojnl (a) . 
Thus, accounting for nonlocality in Lo^i^{a) would be 
particularly important for trapped-particle waves whose 
phase velocity varies with time, modulated waves (Pa- 
per HI) being an example. 

Finally, let us also consider effects that can be driven 
by dissipation, in case when plasma is coUisional. Since u 
changes rapidly with small a in the presence of a phase- 
space clump or a hole, slow decay of a will cause fre- 
quency downshifting or upshifting, correspondingly. Yet, 
since /3 depends on how localized Fb{J) is, another effect 
is anticipated, namely, as follows. Notice that a friction 
drag (say, proportional to the particle velocity) can cause 
condensation of the trapped distribution near the bottom 
of the wave potential trough [46] . Hence, peaking of F( J) 
can occur, and j3 can decrease gradually down to minus 
one. This represents a frequency sweeping mechanism 
additional to those considered in Refs. 



2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



(a) 





1.0 


(b) 


■ 


0.5 






°« 0.0 




^ ___ 






-0.5 




■ 


-1.0 




. 



0.0 0.5 1.0 1.5 2.0 2.5 



0.0 0.5 1.0 1.5 2.0 2.5 



FIG. 4: (a) Normalized action j{r), Eq. ( |A1| . (b) Weight 
function g{r), Eq. (A6l. The dashed lines denote the separa- 
trix r — 1 and also the asymptote of g{r). 

VI. CONCLUSIONS 

In this paper, we expand on the results that were pre- 
viously reported in our brief letter [llj , performing a sys- 
tematic study of kinetic nonlinearities in NDR of an adi- 
abatic sinusoidal Langmuir wave in coUisionless plasma. 
Specifically, a general NDR is derived in a nondiffcrcntial 
form, allowing for an arbitrary distribution of trapped 
electrons. The linear dielectric function is generalized, 
and the nonlinear kinetic frequency shift wnl is found an- 
alytically as a function of the wave amplitude a. Smooth 
distributions yield wnl oc ^/a, as usual. However, beam- 
like distributions of trapped electrons result in different 
power laws, or even a logarithmic nonlinearity, which are 
derived as asymptotic limits of the same dispersion re- 
lation. Such beams are formed whenever the phase ve- 
locity changes, because the trapped distribution is in au- 
toresonance and thus evolves differently from the passing 
distribution. Hence, even adiabatic wnlIq-) is generally 
nonlocal. 



Appendix A: Auxiliary functions 

Here, we summarize explicit expressions for the special 
functions used in this paper 51J. The analysis presented 
in the main text does not rely on these expressions per 
se, yet those are needed to calculate the numerical values 
of j*, ip{j*), and >c. 

Specifically, the normalized action j(r) reads as 



4 r (r-l)K(r) + E(r), r < 1, 



V2E(r-i), 



r > 1 



[Fig. l4la)], where K and E are the complete elliptic inte- 
grals of the first and second kind [52] , 



K(C) = / (1 - C sin^ 6l)-i/2 de, (A2) 

Jo 

E(C)= / (1-Csin2 0)1/2^0, (A3) 

^0 



with the following properties: 



K(0) = E(0) = ^/2, E(l) = l. 



(A4) 



In particular, one thereby gets j(l) = A/n (which can be 
found as well without introducing elliptic integrals) , and 



2 f K(r), r < 1, 
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also yielding the oscillation period, with Eq. (18 1. From 
Eq. (15), one obtains then 



ri-2E(r)/K(r), r < 1, 

'lir) = { , , A6) 

1 2r-l-2rE(r-i)/K(r-i), r>l 



[Fig. |4l^b)]. Further, the function ^(j) is calculated as 
\I/(j) == ip{r{j)), with Eq. ([22]) yielding 
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'^ii~) = ~ / 9ir)j'ir)dr. 
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(A7) 



FIG. 5: (a) Function ^(r), Eq. (A8l. (b) Function q{r 



Eq. ( |A9| ). The dashed lines denote the separatrix r — 1 and 
also the asymptote of q{r). 



Hence, Eqs. (A5) and (A6| result in 



4 r (l-r)K(r) + (2r-l)E(r), r < 1, 

ih(r) = — X < 

Stt 1 ri/2[2(l-r)K(r-i) + (2r-l)E(r-i)], r > 1, 



(A8) 



and V'(l) = 4/(37r) in particular [Fig.lsFa)]. 

Finally, Q{j) is calculated as Q{j) — q{r{j)), with 



q{r) = 1 - 2j{r)'4>{r) 



(A9) 



[Fig.^b)]. Hence, 



X = / [1 - 2j{r)ij{r)] j'{r) dr « 0.544. (AlO) 

Jo 
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